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Abstract. In this Note, we describe the stationary equilibria and the asymptotic be- 
haviour of an heterogeneous logistic reaction-diffusion equation under the influence of 
autonomous or time-periodic forcing terms. We show that the study of the asymptotic 
behaviour in the time-periodic forcing case can be reduced to the autonomous one, the 
last one being described in function of the "size" of the external perturbation. Our re- 
sults can be interpreted in terms of maximal sustainable yields from populations. We 
briefly discuss this last aspect through a numerical computation. 

Resume: Analyse de modeles de dynamique de populations sous l'influence 
de perturbations externes. Cette Note a pour objet l'etude des etats station- 
naires et du comportement asymptotique d'equations de reaction-diffusion avec coeffi- 
cients heterogenes en espace, auxquelles nous ajoutons un terme de perturbation station- 
naire ou periodique en temps. Nos resultats peuvent s'interpreter en termes de recolte 
maximale supportable par une population. Nous soulignons cet aspect a l'aide d'un 
calcul numerique . 

1. Introduction 

The purpose of this Note is to study the following model: 

(1.1) Ut = V • (A(x)Vu) + u((jl(x) - u{x)u) - f(ut,x)p £ (u), (t,x) e E+ x tt. 

The reaction-diffusion models of the type Ut — V • (A(x)Vu) + u(/i(x) — u(x)u) correspond 
to the natural extension of the classical Fisher model [3]. They were first introduced by 
Shigesada et al. [8] for population dynamics. Our aim is to understand the asymptotic 
behaviour of the solutions of such models, when we add a time-periodic forcing term 
f(uit,x). With such additional term, this can interpreted as an harvesting model with 
seasonal harvesting. In real-life context this perturbation term can arise when a quota is 
set on the harvesters. 

We make the following assumptions on the coefficients: the diffusion matrix A(x) is 
assumed to be of class C 1,a (with a > 0) and uniformly elliptic; i.e. there exists r > 
such that A(x) > r/jy for all i6fl. The functions /i and v belong to L°°(Q). Moreover, 
we assume that there exist v_ and V such that < y_ < v(x) < V for all x in Q. The function 
/ is 1-periodic in the first variable and belongs to C°(R x Q), and the function p e defines 
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a "regularized threshold": it is a C 1 (M) nondecreasing function such that p s (s) = for 
all s < and p £ (s) = 1 for all s > e. This threshold guarantees the non-negativity of the 
solutions of ( 11. ip . 

Two kinds of domains Q are considered: either Q = M. N or Q is a smooth bounded 
domain of M . We qualify the first case, f2 = R , as the sp-case and the second one as 
the bounded case. Indeed, in the sp-case, we assume that z/(x) and f(s,x) 

depend on the variables x = (xi, • • • , xjy) in a space-periodic fashion (i.e. for Li, . . . , L N 
fixed positive numbers, a function g is said to be sp-periodic if g(x + k) = g(x) for all 
x G M N and G L\L X • • • x L N Z). In the bounded case, throughout this paper, we 
assume that we have Neumann boundary conditions on dQ. 

2. The case of autonomous forcing 

All the results of this section remain true either in the sp-periodic or bounded cases. 
The proofs are detailed in [B|. 

We consider the equation (11.11) with f{wt,x) = 5h(x), i.e. 

(2.1) u t = V • (A(x)Vu) +u(n(x) -v(x)u) - Sh(x)p e (u), (t,x) e R+ x Q, 

where h is a continuous function such that there exist a,/3 > with a < h(x) < 
(3 for all i6fl, and which is sp-periodic in the sp-case. 

Let Ai be defined as the unique real number such that there exists a function <p > 
which satisfies 

(2.2) - V • (A(x)V(/)) - p{x)(f) = Ai0 in Q, <p > and = 1, 

with either periodic or Neumann boundary conditions, depending on Q, as mentioned 
above. The function <fr is uniquely defined by (12.21) (the existence and uniqueness of Ai 
and follow from the standard Krein-Rutman theory). 

Remark 2.1. Note that if we assume that Ai < and 5 = 0, then, given any continuous 
and bounded function u , the solution u(t, x) of (12. ip with initial data u converges to a 
function p which is the unique bounded and positive solution of V • (A(x)'Vp) +p(p(x) — 
v(x)p) = 0, x G Q. These convergence, as well as existence and uniqueness results are 
proved in [I]. 

We first describe the steady states of (12.11) without "regularized threshold": 

(2.3) V • (A(x)Vps) +ps(p(x) - is(x)ps) - 6h(x) = 0, x G £1. 

Using a Leray-Schauder degree argument, together with the uniqueness of the solution p 
defined in the above remark, we prove the following 

Theorem 2.1. There exists S* > such that for all S s. t. < 5 < S* , (12. 3p admits two 
distinct positive solutions, pj and p\. Moreover, p\ — > and p\ — >■ p uniformly in Q as 
rS - > (I. 

A?0 A? 

Let us set <p := min0(x), 5\ := — — ; — = — — and 62 '■= • Then we have the following 

— xen fju(l + <py Aau_ 

theorem: 
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Figure 1. N — 2, A — I 2 , v = 1, h = 1 and fj,(x) is (1, l)-periodic and 
takes two values, fi = — 1 on f2_ and fi = 10 on where f2 + consists on A; 2 
equally-spaced disks such that, on each period cell [0, l] 2 , |f2 + H[0, 1] 2 | = 1/2, 
and Q~ = M. 2 \Q + . (a) A period cell with k = 2 (b) k = 10; Q + is represented 
in black, (c) The values of Si (continuous line) and S 2 (dashed line) in 
function of k. 

Theorem 2.2. - (i) If Xi < and S < Si, then there exists a positive bounded solution 

Ps of f )2.3p such that ps > — — -. ^ (in particular maxp^ > ). 

u(l + (p) zv 

(ii) If Ai < and 5 > 62, or if Xi > ; there is no positive bounded solution of (12. 3p . 

The proof relies on monotone methods of sub- and super-solutions. For the existence 
result (i), we have computed a sub-solution of the form K<f) with k > 0. The optimal 
value of k, in the sense that it gives the highest value of 5i, is k = —Xi/(p + u</>). We 
have numerically computed the values of 8\ and 5 2 in several particular examples of sp- 
case (see Figured]). The results illustrate the effect of environmental fragmentation on 
the maximum sustainable yield, and show that the interval (^1,^2] on which we have no 
theoretical information can be very narrow (see Figure [l]-(c)). 

Let us turn to the study of the evolution equation (12. ip . We assume that Ai < and e 

is such that Eq := 2— < -; we prove the following theorem: 

<p 2 

Theorem 2.3. Let u(t,x) be the solution of ( 12. ID with initial data u(0,x) = p(x) defined 
in Remark \2.1\ Then u is non-increasing in t and we have the following asymptotic 
behaviour 

(i) if 5 < Si, u(t, x) — > ps(x) uniformly in Q as t +00, where p$ is the unique positive 
maximal solution of f)2.3jl / and 

(ii) if 5 > 5 2 , then u(t,x) < e fort large enough. 

In the above theorem, we assume that u(0, x) = p(x). This means that harvesting starts on 
a stabilized population governed by the standard Fisher model without external forcing. 
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Remark 2.2. These results are sharper than those which could be obtained by a standard 
La Salle invariance principle, since we obtain here discriminatory bounds on 5, which 
determine the asymptotic behaviour of the solutions. 



In this section we consider the general equation ( 11. ip in the bounded case, with 
oj > defined as the frequency of the forcing term. All the results are proved in [2]. 
Let us introduce T := u^ 1 . It is known that under the above assumptions on A(x), 
Au = — V-(A(x)Vu) is a sectorial operator with domain V {A) = {nG H 2 (Q), s. t. d n u = 
on dQ} (see e.g. |5J). As a consequence, —A generates an analytic semigroup e~ At on 
L 2 (fl). Let {^ 2r }{r>o} be the family of interpolation spaces generated by the fractional 
powers of A, where V 2r = V(A r ) (see [7] for details). The existence of a T-periodic 
solution of the equation ( II. ip can be reached by several procedures (e.g. averaging method 
[1]). We present here a result on the existence of a hyperbolic T-periodic solution, which is 
related to the robustness of a hyperbolic stationary solution of the autonomous equation 
( 12. ip . with 8h(x) = f Q f(s,x)ds. More precisely, 

Theorem 3.1. Assume that equation ( 12. ip has a hyperbolic stationary solution q G V 2r , 
< r < 1. Then there exists u* > such that for every uj > u* , the problem ( 11. ip 
possesses a hyperbolic T-periodic solution u u (t) such that for any t G [0, T], u^it) lives in 
a V 2r -neighborhood of q. Furthermore if u — > +oo, then u^t) —>• q in V 2r . 

The proof uses similar arguments as the one of j7] Theorem 76.1, and is therefore based 
on a Lyapunov- Perron type argument; the existence of such a hyperbolic periodic orbit is 
achieved via a fixed point argument on the following operator: 



where L = A — (/i(x) — 2u(x)q)T, E(q,^) = — 2i/(x)£ 2 , £ belongs to a subset of 
L°°(]R, V 2r ) nC°(M, V 2r ), and P and Q are the associated projectors with the exponential 
dichotomy for equation (12.11) related to the existence of a hyperbolic stationary solution. 

The main interest of Theorem 13 .11 is that it gives a simple sufficient condition to ensure 
the existence of a T-periodic solution of (II. ip and that it allows to localize in physical 
space where this solution can appear. Another interesting aspect of this theorem is that 
it reduces the study of existence and stability of a T-periodic solution of (II. ip to that 
of the hyperbolic equilibria of the autonomous version (12.11) . For instance we get as 
an application of Theorems 2 and 4 for Ai < and f(wt,x) = 5g(wt,x) with a < 

g(s, x)ds < (3 for all x G f2, that if 5 < 5\ and u sufficiently large, then there exists a 
stable non-trivial T-periodic solution u^^s of (II. ip in a neighborhood of a solution p§ of 



3. Time-periodic forcing 





(Q. 
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Note added for this version on ArXiV: The reference j6] cited here, corresponds 
the article of the authors entitled "On Population Resilience to External Perturbations" 
which has been published in SI AM J. Appl. math (SIAP), 68 (1), (2007) 133-153. We 
kept here the old reference [6] as in the original article published in C. R. Acad. Sci. 
Paris, Ser. I, 343: 307-310; before the SIAP article. The reference [2] below is still under 
preparation. 
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